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Abstract. In the context of three-neutrino oscillations, we study the possibility of using antineutrinos 
from nuclear reactors to explore the 10“^ eV^ < ArriQ ^8x10 ^ eV^ region of the LMA MSW solu¬ 
tion of the solar neutrino problem and measure Amg with high precision. The KamLAND experiment 
is not expected to determine Amg if the latter happens to lie in the indicated region. By analysing 
both the total event rate suppression and the energy spectrum distortion caused by Pg oscillations in 
vacuum, we show that the optimal baseline of such an experiment is L ~ (20 — 25) km. Furthermore, 
for 10“^ eV^ < Amg ^ 5 x 10“^ eV^, the same experiment might be used to try to distinguish 
between the two possible types of neutrino mass spectrum - with normal or with inverted hierarchy, 
by exploring the effect of interference between the atmospheric- and solar- Am^ driven oscillations; 
for larger values of Amg not exceeding 8.0 x 10“^ eV^, a shorter baseline, L = 10 km, would be 
needed for the purpose. The indicated interference effect modifies in a characteristic way the energy 
spectrum of detected events. Distinguishing between the two types of neutrino mass spectrum requires, 
however, a high precision determination of the atmospheric Am^, a sufficiently large sin^ 9 and a non- 
maximal sin^ 20g, where 9 and 9q are the mixing angles respectively limited by the CHOOZ and Palo 
Verde data and characterizing the solar neutrino oscillations. It also requires a relatively high precision 
measurement of the positron spectrum in the reaction Ug + p ^ e'^ + n. 
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1 Introduction 


In recent years the experiments with solar and atmospheric neutrinos collected strong evidences 
in favor of the existence of oscillations between the flavour neutrinos, and Further progress in 

our understanding of the neutrino mixing and oscillations requires, in particular, precise measurements 
of the parameters entering into the oscillation probabilities - the neutrino mass-squared differences and 
mixing angles, and the reconstruction of the neutrino mass spectrum. 

The atmospheric neutrino data can be explained by dominant ^ v^- and —> Vr oscillations, 
characterized by large, possibly maximal, mixing, and a mass squared difference, having a 

value in the range ||T]] (99% C.L.): 

1.3 X IQ-^eV^ ^ ^ 5 x IQ-^eV^. (1) 


The first results from the Sudbury Neutrino Observatory (SNO) [Q], combined with the mean event rate 
data from the Super-Kamiokande (SK) experiment [|1], provide a very strong evidence for oscillations 
of the solar neutrinos 0 - [IF]. Global analyses of the solar neutrino data, including the SNO results 
and the SK data on the e“—spectrum and day-night asymmetry, show that the data favor the large 
mixing angle (LMA) MSW solution of the solar neutrino problem, with the corresponding neutrino 
mixing parameter sin^ 29q and mass-squared difference Attiq lying in the regions (99.73% C.L.): 


2 X IQ-^eV^ ^ Am| ^ 8 x IQ-'^eV^ 
0.6 < sin^ 200 < 1. 


( 2 ) 

( 3 ) 


The best fit value of Attiq found in the independent analyses ^ is spread in the interval 

(4.3 — 6.3) X 10“^ eV^. The results obtained in [^, ^ show that values of Amg > 10“^ eV^ are 

allowed already at 90% C.L. Values of cos 20 q < 0 (for Amg > 0) are disfavored by the data. 

Important constraints on the oscillations of electron (anti-)neutrinos, which play a significant role 
in our current understanding of the possible patterns of oscillations of the three flavour neutrinos and 
anti-neutrinos, were obtained in the CHOOZ and Palo Verde disappearance experiments with reactor tJ^ 
00]- The CHOOZ and Palo Verde experiments were sensitive to values of neutrino mass squared 
difference Am^ > 10“^ eV^, which includes the corresponding atmospheric neutrino region, eq. (1). 
No disappearance of the reactor Pg was observed. Performing a two-neutrino oscillation analysis, the 
following rather stringent upper bound on the value of the corresponding mixing angle, 9, was obtained 
by the CHOOZ collaboration [] | [TT[ ] at 95% C.L. for Am? > 1.5 x 10“^eV^: 

sin^ 9 < 0.09. (4) 


The precise upper limit in eq. (^ is Am^-dependent: it is a decreasing function of Am? as Am? 
increases up to Am? ~ 6 ■ 10“^ eV^ with a minimum value sin^ 9 ~ 10“^. The upper limit becomes 
an increasing function of Am^ when the latter increases further up to Am? ~ 8 • 10“^ eV^, where 
sin^ 0 < 2 • 10“^. Somewhat weaker constraints on sin^ 0 have been obtained by the Palo Verde 
collaboration In the future, sin^ 0 might be further constrained or determined, e.g., in long baseline 


neutrino oscillation experiments []13|]. 


The long baseline experiment with reactor Pg Kam i. AND [ |I^ has been designed to test the LMA 
MSW solution of the solar neutrino problem. This experiment is planned to provide a rather precise 
measurement of Amg and sin^ 200. Due to the long baseline of the experiment, L ~ 180 km, however, 
Amg can be determined with a relatively good precision only if Amg ^ 10“^ eV^. 


*The possibility of large sin^ 0 > 0.9 which is admitted hy the CHOOZ data alone is incompatihle with the neutrino 
oscillation interpretation of the solar neutrino deficit (see, e.g., [p^ p^]) 
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The explanation of both the atmospheric and solar neutrino data in terms of neutrino oscillations 
requires, as is well-known, the existence of 3-neutrino mixing in the weak charged lepton current: 


t'iL = ^ UijUji^, 


( 5 ) 


i=i 


where un, I = e, r, are the three left-handed flavour neutrino fields, is the left-handed field of 
the neutrino Vj having a mass rrij > 0 and [/ is a 3 x 3 unitary mixing matrix - the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) neutrino mixing matrix [ |l7| , |T^. The three neutrino masses mi_ 2,3 can 
obey the so-called normal hierarchy (NH) relation mi < m 2 < m3, or that of the inverted hierarchy 
(IH) type, m 3 < mi < m 2 . Thus, in order to reconstruct the neutrino mass spectrum in the case 
of 3-neutrino mixing, it is necessary to establish, in particular, which of the two possible types of 
neutrino mass spectrum is actually realized. This information is particularly important for the studies 
of a number of fundamental issues related to lepton mixing, as like the possible Majorana nature of 
massive neutrinos, which can manifest itself in the existence of neutrino-less double /3-decay (see, e.g., 
[jT^ 0]). It would also constitute a critical test for theoretical models of fermionic mass matrices and 


flavor physics in general. 

It would be possible to determine whether the neutrino mass spectrum is with normal or inverted 
hierarchy in terrestrial neutrino oscillation experiments with a sufficiently long baseline, so that the 
neutrino oscillations take place in the Earth and the Earth matter effects in the oscillations are non- 


negligible |21, 22, 23]. The ambiguity regarding the type of the neutrino mass spectrum might be 


resolved by the MINOS experiment [|13|], although on the baseline of this experiment the matter effects 


are relatively small [21]. This might be done in an experiment with atmospheric neutrinos, utilizing a 
detector with a sufficiently good muon charge discrimination [^]. The experiments at neutrino factories 


would be particularly suitable for the indicated purpose [22, 23]. 

In this paper, in the context of three-neutrino oscillations, we study the possibility of using anti¬ 
neutrinos from nuclear reactors to explore the Amg > 10“^ eV^ region of the EMA MSW solution. 
Such an experiment might be of considerable interest if, in particular, the results of the KamEAND 
experiment will confirm the validity of the EMA-MSW solution of the solar neutrino problem, but will 
allow to obtain only a lower bound on Amg due to the fact that Amg > 10 “^ eV^ [^, 

We determine the optimal baseline of the possible experiment with reactor which would provide a 
precise measurement of Amg in the region 10“^ eV^ < Amg ^ 8 x 10“^ eV^. Eurthermore, the 
same experiment might be used to try to distinguish between the two types of neutrino mass spectrum 
- with normal or with inverted hierarchy. This might be done by exploring the effect of interference 
between the amplitudes of neutrino oscillations, driven by the solar and atmospheric Am^, i.e., by 
Amg and Eor the optimal baseline found earlier, L = (20 — 25) km, the indicated effect 

could be relevant for 10“^ eV^ < Amg ^ 5 x 10“^ eV^. Eor larger values of Arrig 
interval 


< Amg ^ 5 X 10““* eV^. Eor larger values of Arrig within the 
), the effect could be relevant at L = 10 km. Distinguishing between the two possible types 
of neutrino mass spectrum requires a relatively high precision measurement of the positron spectrum 
in the reaction + p —i e"'' -|- n (i.e., a high statistics experiment with sufficiently good energy 
resolution), a measurement of Am^^^ with very high precision, sin^ 20g 7 ^ 1.0, e.g., sin^ 20g ^ 0.9, 
and a sufficiently large value of the angle 9, which for Amg <C Am^^^ controls, e.g., the oscillations 
of the atmospheric v^. and and is constrained by the CHOOZ and Palo Verde data. 

2 The zzg Survival Probability 


We shall assume in what follows that the 3-neutrino mixing described by eq. (^) takes place. We 
shall number (without loss of generality) the neutrinos with definite mass in vacuum Vj, j = 1, 2, 3, 
in such a way that their masses obey mi < m2 < m3. Then the cases of NH and IH neutrino mass 
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spectrum differ, in particular, by the relation between the mixing matrix elements \Uej\, j = 1,2, 3, and 
the mixing angles 6 q and 9 (see further). With the indicated choice one has > 0 for j > k. Let 

us emphasize that we do not assume any of the relations mi <C m 2 'C m 3 , or mi ^ m 2 4; m 3 , or 
mi <C m 2 = m 3 , to be valid in what follows. 

Under the conditions of the experiment we are going to discuss, which must have a baseline L 
considerably shorter than the baseline ~ 180 km of the Kam i. AND experiment, the reactor Ug oscilla¬ 
tions will not be affected by Earth matter effects when the travel between the source (reactor) and 
the detector. If 3-neutrino mixing takes place, eq. (^), the Ug would take part in 3-neutrino oscillations 
in vacuum on the way to the detector. 

We shall obtain next the expressions for the reactor Pg survival probability of interest in terms 
of measurable quantities for the two types of neutrino mass spectrum. In the case of normal hierarchy 
between the neutrino masses we have: 


Am| = Amh, 


and 


|Uel| = cos 001/1 - |(7e3p, \Ue2\ = sin6»0y/l - |C/e3p, 

9q = 612 , |Ue 3 p = sin^ 9 = sin^ 9i 


( 6 ) 

(7) 


where 

" hs, ( 8 ) 

012 and 013 being two of the three mixing angles in the standard parameterization of the PMNS matrix 
(see, e.g., [1^1). Note that |C 7 e 3 |^ is constrained by the CHOOZ and Palo Verde results. It is not difficult 
to derive the expression for the i?g survival probability in the case under discussion: 

PNH{i>e ^e) 

= 1 — 2 sin^ 0 cos^ 0(1 — cos ^ 


— - cos^ 0 sin^ 200 ( 1 — cos 


2 E^ 
AmQ L 
2 E^ 


+ 2 sin^ 0 cos^ 0 sim 00 ( cos , , 

' ' 2 E^, 2E^ J 2E„ 

where E^^ is the neutrino energy and we have made use of eqs. (|), ( 0 ) and (|). 


Amoi L Am%L\ Amo, L 


If the neutrino mass spectrum is with inverted hierarchy one has (see, e.g., ||^, ^ |160): 

Am% = Am| 2 , 


and 


\Ug2\ = COS 00 V 1 - |i7el|2, |(7e3| = sin 0q ^/l - |C/el P. 

The mixing matrix element constrained by the CHOOZ and Palo Verde data is now |?7ei : 

|Uei|^ = sin^ 0 . 

The expression for the Ug survival probability can be written in the form : 

Piuii'e Pe) 

= l-2sin2 0cos2 0 ri-cos-^""''^ 


— - cos^ 0 sin^ 200 ( 1 — cos 


2 E^ 
AmQ L 
2 E„ 


+ 2 sin^ 0 cos^ 0 cos^ 00 ( cos 


I cos I - 


/ Amqi L Am% L 


^31 

2 E,. 


2 E,. 


— cos 


Am^i L 

2 E„ 


(9) 


( 10 ) 

( 11 ) 

( 12 ) 


( 13 ) 
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Several comments concerning the expressions for the v^. survival probability, eqs. (^) and (p3[), 
follow. In the first lines in the right-hand side of eqs. @ and (|T^), the oscillations of the electron 
(anti-)neutrino driven by the “atmospheric” are accounted for. The CHOOZ and Palo Verde 

experiments are primarily sensitive to this term and their results limit sin^ Q. The second lines in the 
expressions in eqs. @ and ( [T^ contain the solar neutrino oscillation parameters. This is the term 
KamLAND should be most sensitive to. For Amg <C Am^]^ = Amg ^ 10“^ eV^, only one 

of the indicated two terms leads to an oscillatory dependence of the Fg survival probability for the ranges 
oiL/Ey characterizing the CHOOZ and Palo Verde, and the KamLAND experiments: on the source- 
detector distance L of the CHOOZ and Palo Verde experiments the oscillations due to Amg cannot 
develop, while on the distance(s) traveled by the Pg in the KamLAND experiment Am^j^ causes fast 
oscillations which average out and are not predicted to lead, e.g., to specific spectrum distortions of the 
KamLAND event rate. 

The terms in the third lines in eqs. and ( p^ are not present in any two-neutrino oscillation 
analysis. They represent interference terms between the amplitudes of neutrino oscillations, driven 
by the solar and atmospheric neutrino mass squared differences. The term in eq. is proportional 
to sin^ 0g, while the corresponding term in eq. ( pA| ) is proportional to cos^ 0g [|2^. This is the only 
difference between P/v/t(Fg —> Pg) and Pinii^e Pe)^ that can be used to distinguish between the 
two cases of neutrino mass spectrum in an experiment with reactor Pg. Obviously, if cos 20g = 0, we 
have PNuii'e —> Fg) = Pini^e Pe) and the two types of spectrum would be indistinguishable in the 
experiments under discussion. For vanishing sin^ 9, only the terms in the second line of eqs. (^) and 
© survive, and the two-neutrino mixing formula for solar neutrino oscillations in vacuum is exactly 
reproduced. 

Let us discuss next the ranges of values the different oscillation parameters, which enter into 
the expressions for the probabilities of interest P^nii^e Pe) and Pjnii^e —> Pe); can take. The 
allowed region of values of Amg^, Amg, sin^ 6 q and 6 should be determined in a global 3-neutrino 
oscillation analysis of the solar, atmospheric and reactor neutrino oscillation data, in which, in par¬ 
ticular, Amg should be allowed to take values in the LMA solution region, including the interval 
Amg ~ (1.0 — 6.0) X lO”"^ eV^. Such an analysis is lacking in the literature. However, as was shown 
in [0Ol], a global analysis of the indicated type would not change essentially the results for the LMA 
MSW solution we have quoted | in eqs. and (|3|) as long as Am|]^ 1.5 x 10“^ eV^. The reason is 

that for Amg^ ^ 1.5 x 10“^ eV^ and Amg 6.0 x 10“^ eV^, the solar survival probability, which 
determines the level of suppression of the solar neutrino flux and plays a major role in the analyses of 
the solar neutrino data, depends very weakly on (i.e., is practically independent of) Am|]^. Thus, Amg 
and 9q are uniquely determined by the solar neutrino and CHOOZ and Palo Verde data, independently 
of the atmospheric neutrino data and of the type of the neutrino mass spectrum. The CHOOZ and 
Palo Verde data lead to an upper limit on Amg in the LMA MSW solution region (see, e.g., [^, |^): 
Amg ^ 7.5 X 10“^ eV^. For Amg ^ 1.0 x 10“^ eV^, the CHOOZ and solar neutrino data imply the 
upper limit on sin^ 6 given in eq. (^. For Amg ~ (2.0 — 6.0) x 10“^ eV^ of interest, the upper limit 
on sin^ 0 as a function of Am|^ ^ 10“^ eV^ for given Amg and sin^ 20g is somewhat more stringent 

iHi- 

Would a global 3-neutrino oscillation analysis of the solar, atmospheric and reactor neutrino 
oscillation data lead to drastically different results for Am|g in the two cases of normal and inverted 
neutrino mass hierarchy? Our preliminary analysis shows that given the existing atmospheric neutrino 
data from the Super-Kamiokande experiment, such an analysis i) would not be able to discriminate 
between the two cases of neutrino mass spectrum, and ii) would give essentially the same allowed 
region for Amgj^ in the two cases of neutrino mass spectrum. We expect the regions of allowed values 

^Let us note that the LMA MSW solution values of ArriQ and Oq we quote in eqs. (^) and were obtained hy taking 
into account the CHOOZ and Palo Verde limits as well. 


4 




of the mixing angle Oatm^ which controls the dominant atmospheric and ^ Vj. oscillations, 

to differ somewhat in the two cases. Note, however, that this mixing angle does not enter the expression 
for the Vf. survival probability we are interested in. 


For Attiq ^ 1.0 X 10 ^ eV^ and sufficiently small values of sin^ 


coincides effectively 

of the two-neutrino and oscillation analyses of the SK atmospheric neutrino data. 
If sin^ 6 > 0.01, a three-neutrino oscillation analysis of the atmospheric neutrino and CHOOZ data. 


with Am^ 


performed under the assumption of Attiq 1.0 x 10“^ eV^ [^], gives regions of allowed values of 
Am^^rrx = Aml]^, which are correlated with the value of sin^ 6 . The latter must satisfy the CHOOZ 
and Palo Verde constraints. 


At present, as we have already indicated, a complete three-neutrino oscillation analysis of the 
atmospheric neutrino and CHOOZ data with Attiq allowed to take values up to ~ (6.0 — 7.0) x 
10“^ eV^, i.e., in the region where deviations from the two-neutrino approximation could be non- 
negligible, is lacking in the literature. Therefore in what follows we will use representative values of 
Aml^ which lie in the region given by eq. (|]). 


3 The Difference between PNH{j^e and ^e) 


Let us discuss next in greater detail the difference between the Fg surviving probabilities in the 
two cases of neutrino mass spectrum of interest, PNHp'e Pe) and Pinii^e Pe)- While the terms 
in the first two lines in eqs. and ( p^ describe oscillations in LjEy with frequencies and 

Amg/dvr, respectively, the third term has the shape of beats, being produced by the interference of two 
waves, with the same amplitude but slightly different frequencies: 


cos 


/ Amj^L 

V 


L\ Arriqi L 

- - cos - 

2 Ey J 2 Ey 




ArrP L 

2 sin -—— sin 

4.Ey 

Amg L 

2 sin -— sin 

4.Ey 


L 

2Ey “ 

Am|^\ 
2Ey J 


AttIq 

4K J 


(14) 


This is a modulated oscillation with approximately the same frequency of the first term in eqs. @ and 
( p^ (Arngj^/dvr) and amplitude oscillating between 0 and 2sin^ 6 q of the amplitude of the first term 
itself. The beat frequency is equal to the frequency of the dominant oscillation (Amg/dvr). The modu¬ 
lation is exactly in phase with the Amg—driven dominant oscillation of interest, so that the maximum of 
the oscillation amplitude of the interference term (third lines in the expressions for PnhP'e ^e) and 
Pini^e Fe)) is reached in coincidence with the points of maximal decreasing of the Fg survival prob¬ 
ability, where Amg L/AE = 7r/2, and vice versa - this amplitude vanishes at the local maxima of the 
survival probability. At the minima of the Fg survival probability, for instance at Amg L/AEy = 7r/2, 
pNH(iH) takes the value: 






= 1 — 2 sin^ 9 cos^ 9 — cos^ 9 sin^ 29p 


=1 


© 


cos 29Q 2 sin^ 9 cos^ 9 cos . 


(15) 


From eqs. (^), ( pA| ) and one deduces that: 

• for maximal mixing, cos 29q = 0, the last term cancels, and 

• for very small mixing angles, cos 20g ~ 1, the terms describing the oscillations driven by Am|g 
in the NH and IH cases have opposite signs: the two waves are exactly out of phase. 
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for intermediate values of cos 200 from the LMA MSW solution region, cos 29q, = (0.3 — 0.6), 
the driven contributions in the cases of normal and inverted hierarchy have still opposite 

signs and the magnitude of the effect is proportional to 2 cos 29q sin^ 9. 


The net result of these properties is that in the region of the minima of the survival probability 
due to Attiq, where Amg L / {2E) = 7r(2/c + 1), A: = 0,1,..., the difference between ^ i^e) 

and Pjnii'e —> ^e)) is maximal. In contrast, at the maxima of P^nii^e i^e) and Piuii^e —^ f'e)) 
determined by Am-Q L / {2E) = 2TTk, we have, for any sin^ 9q, ENui^'e —> i^e) = Piuik'e i^e)- 

The two-neutrino oscillation approximation used in the analysis of the CHOOZ and Palo Verde 
data is rather accurate as long as AttIq is sufficiently small [p^: for AttIq ^ 10“^ eV^, the L/Ey 
values characterizing these experiments, chosen to ensure maximal sensitivity to Am|^ ^ 10“^ eV^, 
are much smaller than the value at which the first minimum of PNH{iH){^e Pe) due to the AtUq- 
dependent oscillating term occurs. Correspondingly, the effect of the interference term is strongly 
suppressed by the beats. For Auiq ^ 2 x 10“^ eV^ this is no longer valid and the interference term 


under discussion has to be taken into account in the analyses of the CHOOZ and Palo Verde data |29]. 


4 Measuring Large Am| at Reactor Facilities 


As is well-known, nuclear reactors are intense sources of low energy Pg {Ey ^ 8 MeV), emitted 
isotropically in the /3-decays of fission products with high neutron density [^]. Anti-neutrinos can then 
be detected through the positrons produced by inverse /3-decay on nucleons. The reactor Pg energy spec¬ 
trum has been accurately measured and is theoretically well understood ^ : it essentially consists 

of a bell-shaped distribution in energy centered around Ey ^ 4: MeV, having a width of approximately 
3 MeV. CHOOZ, Palo Verde and KamLAND are examples of experiments with reactor Pg, the main 
difference being the distance between the source and the detector explored (L ~ 1 km for CHOOZ and 
Palo Verde, and L ~ 180 km for KamLAND). 


The best sensitivity to a given value of Ami, of the experiment of interest is at L at which the 
maximum reduction of the survival probability is realized. As can be seen from eqs. (0) - ([T3|), this 
happens for L around L* = 2tt Ey/ Auiq. This implies that for Ey = 4 MeV, the optimal length to test 
neutrino oscillations with reactor experiments is: 


L* ^ 


5 X IQ-^ 
(Am|/eV^) 


(16) 


The best sensitivity of KamLAND, for instance, is in the range of 2 3 x 10“® eV^. We will discuss 

next in greater detail the distances L which could be used to probe the LMA MSW solution region at 
Am| > 10“^ eV^, in order to extract Am| from these oscillation experiments. 

4.1 Total Event Rate Analysis 

One of the signatures of the Pg—oscillations would be a substantial reduction of the measured to¬ 
tal event rate due to the reactor Pg in comparison with the predicted one in the absence of oscillations. In 
order to compute the expected total event rate one has to integrate the Pg survival probability multiplied 
by the Pg energy spectrum over Ey. In Fig. we show this averaged survival probability for different 
values of L as a function of Am|, using the “best fit” values ^ ^ ^ for Am|^ and sin^ 200. When 
averaging over the Pg energy spectrum, oscillatory effects with too short a period are washed out, and 

^By reactor Ue energy spectrum we mean here and in what follows the product of the i/e production spectrum and the 
inverse /3-decay cross-section, which gives the “detected” neutrino spectrum in the no oscillation case. The i/e production 
spectrum is known with larger uncertainties at tie energies Ey ^ 2 MeV, but this range is not of interest due to the threshold 
energy = 1.8 MeV of the inverse /3-decay reaction [p^]. Certain known time dependence at the level of a few percent is 
also present up to 3.5 MeV [p5|| and should possibly be taken into account in the analysis of the experimental data. 
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0.1 0.2 0.5 1 2 5 10 

Am^ (10”^ eV^) 


Figure 1: The reactor Fg survival probability, averaged over the Fg energy spectrum, for = 

2.5 X 10“^ eV^, sin^ 20q = 0.8, sin^ 6 = 0.05, as a function of Atuq. The curves correspond to 
L = 180 km (long dashed), L = 50 km (dashed), L = 20 km (thick) and L = 10 km (dotted), 
respectively. 


the experiment is sensitive only to the average amplitude. This happens when the width 5Ey of the 
energy spectrum is such that the integration runs over more than one period, i.e., approximately for: 


^ All El dxlO^eV^ 

A ^ _ 

^ Am?'L Am^{L/Km) 


(17) 


Since 6Eiy ~ 3 MeV, at KamLAND this happens approximately for Attiq ^ 7 x 10“® eV^. The 
corresponding curve in Fig. 1 indicates that the actual sensitivity extends to somewhat larger values 
of Attiq than what is expected on the basis on the above estimate, but the total event rate becomes 
flat for Atttq ^ 10“'^ eV^. This means that KamLAND will be able, through the measurement of the 
total even rate, to test all the region of the LMA MSW solution and determine whether the latter is the 
correct solution of the solar neutrino problem, but will provide a precise measurement of Amg only if 
Attiq ^ 10“^^ eV^. If AmQ ^ 2 x 10“^ eV^, it would be possible to obtain only a lower bound on 
Arrig and a new experiment might be required to determine Amg. 

Fig. 1 shows that as L decreases, the sensitivity region moves to larger Amg. These results imply that 
a reactor Fg experiment with L = (20 — 25) km can probe the range 0.8 x 10“^ eV^ < Amg ^ 6 x 
10“^ eV^. One finds that for Amg = 2 x lO”"^ eV^ and Am'^^ = 2.5 x 10“^ eV^, the best sensitivity is 
at L = 20 km. Moreover, with L = (20 — 25) km, the predicted total event rate deviates from being flat 
(in Amg) actually for Amg as large as ~ (5 — 6) x 10“^ eV^. In order to have a precise determination 
of Amg with L = (20 — 25) km for the largest values given in eq. (^, Amg = (7-^ 8) x lO”"^ eV^, one 
should use the information about the e+—spectrum distortion due to the Fg—oscillations. By measuring 
the e+—spectrum with a sufficient precision it would be possible to cover the whole interval 


1.0 X lO""^ eV^ ^ Am% ^ 8.0 x lO""^ eV^ , (18) 

i.e., to determine Amg if it lies in this interval, by performing an experiment at L = (20 — 25) km 
from the reactor(s) ^ (see the next sub-section). 

''The fact that if Ami, = 3.2 x 10“'' eV^, a reactor Ue experiment with L = 20 km would allow to measure Ami with 
a high precision was also noticed recently in [§]• 
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Figure 2: The reaetor Fe energy speetrum at distance L = 20 km from the source, in the absence of 
Fg oscillations (double-thick solid line) and in the case of Fg oscillations characterized by = 

2.5 X 10“^ eV^, sin^ 26q = 0.8 and sin^ 6 = 0.05. The thick lines are obtained for Aitiq = 2 x 10“^ 
eV^ and correspond to NH (light grey) and IH (dark grey) neutrino mass spectrum. Shown is also the 
spectrum for Ami, = 6 x 10“^ eV^ in the NH (dotted) and IH (dashed) cases. 


Applying eq. (17) with Am^ = Am|^, one sees that for the ranges of L which allow to probe 
Am| from the LMA MSW solution region, the total event rate is not sensitive to the oscillations driven 
by Arri^i ^ 1.5 x 10“^ eV^. Thus, the total event rate analysis would determine Ami which would 
be the same for both the normal and inverted hierarchy neutrino mass spectrum. 

4.2 Energy Spectrum Distortions 


An unambiguous evidence of neutrino oscillations would be the characteristic distortion of the 
Fg energy spectrum. This is caused by the fact that, at fixed L, neutrinos with different energies reach 
the detector in a different oscillation phase, so that some parts of the spectrum would be suppressed 
more strongly by the oscillations than other parts. The search for distortions of the Fg energy spectrum 
is essentially a direct test of the Fg oscillations. It is more effective than the total rate analysis since it 
is not affected, e.g., by the overall normalization of the reactor Fg flux. However, such a test requires a 
sufficiently high statistics and sufficiently good energy resolution of the detector used. 

Energy spectrum distortions can be studied, in principle, in an experiment with L = (20 — 25) 
km. In Fig. ^ we show the comparison between the Fg spectrum expected for Am| = 2 x 10“^ eV^ 
and Am| = 6 x 10“^ eV^ and the spectrum in the absence of Fg oscillations. No averaging has been 
performed and the possible detector resolution is not taken into account. The curves show the product 
of the probabilities given by eqs. (§) and (|I^) and the predicted reactor Fg spectrum [^]. As Fig. 

^ illustrates, the Fg spectrum in the case of oscillation is well distinguishable from that in the absence 
of oscillations. Moreover, for Am| lying in the interval 10“^ eV^ < Am| ^ 8.0 X 10 ^ eV^, the 
shape of the spectrum exhibits a very strong dependence on the value of Am|. A likelihood analysis 
of the data would be able to determine the value of Am| from the indicated interval with a rather good 
precision. This would require a precision in the measurement of the e"''—spectrum, which should be 
just not worse than the precision achieved in the CHOOZ experiment and that planned to be reached in 




the KamLAND experiment. If the energy bins used in the measurement of the spectrum are sufficiently 
large, the value of Atoq thus determined should coincide with value obtained from the analysis of the 
total event rate and should be independent of 

5 Normal vs. Inverted Hierarchy 


In Fig. ^ we show the deformation of the reactor Fg spectrum both for the normal and inverted 
hierarchy neutrino mass spectrum: as long as no integration over the energy is performed, the deforma¬ 
tions in the two cases of neutrino mass spectrum can be considerable, and the sub-leading oscillatory 
effects driven by the atmospheric mass squared difference (see the first and the third line of eqs. - 
dlSl)) can, in principle, be observed. They could be used to distinguish between the two hierarchical pat¬ 
terns, provided the solar mixing is not maximal sin^ 9 is not too small and Am|]^ is known with high 
precision. It should be clear that the possibility we will be discussing next poses remarkable challenges. 


The experiment under discussion could be in principle an alternative to the measurement of 
the sign of in long (very long) baseline neutrino oscillation experiments [^, or in the 

experiments with atmospheric neutrinos (see, e.g., [p^). 

The magnitude of the effect of interest depends, in particular, on three factors, as we have already 
pointed out: 


• the value of the solar mixing angle 6q: the different behavior of the two survival probabilities 
is due to the difference between sin^ 0 q and cos^ 9 q; correspondingly, the effect vanishes for 
maximal mixing; thus, the more the mixing deviates from the maximal the larger the effect; 

• the value of sin^ 9, which controls the magnitude of the sub-leading effects due to Am|]^ on the 
Attiq— driven oscillations: the effect of interest vanishes in the decoupling limit of sin^ 0 —> 0; 

• the value of Ami, (see Fig. |T]): for given L and Am| the difference between the spectrum in the 
cases of normal and inverted hierarchy is maximal at the minima of the survival probability, and 
vanishes at the maxima. 


A rough estimate of the possible difference between the predictions of the event rate spectrum 
for the two hierarchical patterns, is provided by the ratio between the difference and the sum of the two 
corresponding probabilities at Am|L = 2ttE^: 

Pnh — PiH 2 cos 200 sin^ 0 cos^ 0 Am| 

_ — _ COS TT _(19) 

Pnh + PiH 1 — 2 sin^ 0 cos^ 0 — cos^ 0 sin^ 20q Am| 

The ratio could be rather large: the factor in front of the cos tt Am 3 ^/Am| is about 25% for sin^ 200 = 
0.8 and sin^ 0 = 0.05. 

The actual feasibility of the study under discussion depends crucially on the integration over 
(i.e., the binning in) the energy: for the effect not to be strongly suppressed, the energy resolution of 
the detector must satisfy: 

, ^ ^ dvr 2 h- 6 X 10^ eV3 

_ ^ Am j, L ~ Ami, (L/km) ' 

^It would be impossible to distinguish between the normal and inverted hierarchy neutrino mass spectrum if for given 
Atuq > 10““^ eV^ and sin^ 29q ^ 1, the LMA solution region is symmetric with respect to the change 9q 7r/2 - 9q 
(cos 29Q — cos 2 ^ 0 ). While the value of sin^ 29q is expected to he measured with a relatively high precision by the 
KamLAND experiment, the sign of cos 29q will not be fixed by this experiment. However, the 9q — (ti/2 — 9q) ambiguity 
can be resolved by the solar neutrino data. Note also that the current solar neutrino data disfavor values of cos 29q < 0 in the 
LMA solution region (see, e.g., |^, |I^). 


9 













For L ~ 1 km this condition could be satisfied for 6 E 1 , ~ /S.Ey, but at L = (15 — 20) Km, for 
= 2.5 X 10“^ eV^ and Ey in the interval (3 — 5) MeV, one should have /S.Ey ^ 0.5 MeV. 

Our discussion so far was performed for simplicity in terms of the reactor Fg energy spectrum, 
while in the experiments of interest one measures the energy of the positron emitted in the inverse /3- 
decay, Ee- The relation between Ee and Ey is well known (see for instance [36]), and, up to corrections 
of at most few per cent, consists just in a shift due to the threshold energy of the process: Ey = 



Figure 3: Comparison between the predicted event rate spectrum at L = 20 km, measured in energy 
bins having a width of /S.Ey = 0.3 MeV in the cases of normal (light grey) and inverted (dark grey) 
neutrino mass hierarchy. The two upper and the lower left figures are for Atoq = 2 x lO”"^ eV^, 
sin^ 200 = 0.8, sin^ 0 = 0.05, and Am|]^ = 1.3; 2.5; 3.5 x 10“^ eV^, respecfively. The lower righf 
figure was obfained for Attiq = 6 x lO”"^ eV^ and = 2.5 x 10“^ eV^. 

Eg + {El^ — me)- The maximal /S.Ey allowed in order fo make fhe effecl observable can be fhen 
direcfly compared fo fhe experimenfal posifron energy resolufion A£^e 0- 

For Ami, ^ lO”"^ eV^, fhe firsf (mosf significant) minimum of the survival probability can be 
explored if L ~ 180 km. In this case, due to the bigger distance L, the energy resolution required would 
be by a factor of ten smaller. This means that for Am| <C Amg^, it is practically impossible to realize 
the condition of maximization of the difference between the survival probabilities in the two cases of 
neutrino mass spectrum without strongly suppressing the magnitude of the difference by the binning of 
the energy spectrum. 

In order to illustrate what are the concrete possibilities in the case of the experiment under dis¬ 
cussion, we have divided the energy interval 2.7 MeV < Ey < 7.2 MeV into 15 bins, with AEy = 0.3 
MeV, and calculated the value of the product of the survival probability and the energy spectrum in each 
of the bins. The results are shown in Fig. |[ 

®In the CHOOZ experiment, for instance, the binning in Ee was AEe ~ 0.40 MeV Jul]. KamLAND is expected to have 
a resolution better than AEe/Ee = 10%/\/Ee, where Ee is in MeV H- 


10 





























( 21 ) 


As our results show and Fig. 3 indicates, for Arrig^ = (1.5 — 3.0) x 10 ^ eV^, 

Am| ^ (2.0 - 5.0) X lO""^ eV^ 

sin^ 200 = 0.8 and sin^ 9 = (0.02 — 0.05), it might be possible to distinguish the two cases of neutrino 
mass spectrum by a high precision measurement of the positron energy spectrum in an experiment with 
reactor Fg with a baseline of L = (20 — 25) km. This should be a high statistics experiment (not less 
than about 2000 Fe—induced events per year) with a sufficiently good energy resolution Q For larger 
values of Atoq not exceeding 8.0 x 10“^ eV^, and Am|^ = (1.5 — 3.0) x 10“^ eV^, the experiment 
should be done with a smaller baseline, L = 10 km. If, however, sin^ 9 ^ 0.01, and/or sin^ 20q ^ 0.9, 
and/or sin^ 200 ^ 0.9 but the LMA solution admits equally positive and negative values of cos 200 , the 
difference between the spectra in the two cases becomes hardly observable. Further, in obtaining Fig. 3 
we have implicitly assumed that Am^i is known with negligible uncertainty. Actually, for the difference 
between the spectra under discussion to be observable, has to be determined, according to our 

estimates, with a precision of ~ 10 % or better given the values of Atuq, sin^ 200 and sin^ 0 , a 
spectrum in the NH case corresponding to a given Arn^i can be rather close in shape to the spectrum 
in the IH case for a different value of Am|^. There is no similar effect when varying Attiq. 

6 Conclusions 

Reactor experiments have the possibility to test the LMA MSW solution of the solar neutrino 
problem. While the Kam i. AND experiment should be able to test this solution, a new experiment 
with a shorter baseline might be required to determine AtUq with high precision if the results of the 
KamLAND experiment show that Attiq > 10“^ eV^. Performing a three-neutrino oscillation anal¬ 
ysis of both the total event rate suppression and the e+—energy spectrum distortion caused by the 
Fg—oscillations in vacuum, we show that a value of Attiq from the interval 10 “^ eV^ < Attiq ^ 8.0 x 
10“^ eV^ could be determined with a high precision in experiments with L = (20 — 25) km if 
the e"''—energy spectrum is measured with a sufficiently good accuracy. Furthermore, if Aitiq = 
(1.0 — 5.0) X 10“^ eV^, such an experiment with L = (20 — 25) km might also be able to distinguish 
between the cases of neutrino mass spectrum with normal and inverted hierarchy; for larger values of 
Attiq not exceeding 8.0 x 10 “^ eV^, a shorter baseline, L = 10 km, should be used for the pur¬ 
pose. The indicated possibility poses remarkable challenges and might be realized for a limited range 
of values of the relevant parameters. The corresponding detector must have a good energy resolution 
(allowing a binning in the positron energy with AE^ ^ 0.40 MeV) and the observed event rate due to 
the reactor Fg must be sufficiently high to permit a high precision measurement of the e+—spectrum. 
Further, the mixing angle constrained by the CHOOZ and Palo Verde data 0 must be sufficiently large 
(sin^ 0 ~ 0.03 — 0.05), and the “solar” mixing angle 00 should not be maximal (sin^ 200 ^ 0.9). In 
addition, the value of Am|^, which is responsible for the dominant ^ and F^ —> Ur oscillations 
of the atmospheric neutrinos, should be known with a high precision. However, as it is well known, 
“only those who wager can win” 
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’Preliminary estimates show that a detector of the type of KamLAND and a system of nuclear reactors with a total power 
of approximately 5-6 GW might produce the required statistics and precision in the measurement of the positron spectrum. 

*The analysis, e.g, of the MINOS potential for a high precision determination of Am|i in the case of Atuq ^ 10“"* * eV^ 
yields very encouraging results (see, e.g., |^8|]). For Amf = (2.0 — 8.0) x 10“^ eV^, such analysis is lacking in the 
literature. 


11 




References 

[1] Super-Kamiokande Collaboration, Y. Fukuda et al., Phys. Rev. Lett. 81 (1998) 1562. 


[2] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87 (2001) 071301 (nucl-ex/0106015). 

[3] S. Fukuda et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 86 (2001) 5651 (hep- 
ex/01 03032). 


[4] V Barger, D. Marfatia and K. Whisnant, hep-ph/0106207 


[5] G. L. Fogli et al., \hep-ph/010624 7 . 

[6] J. N. Bahcall, M. C. Gonzalez-Garcia and C. Pena-Garay, JHEP 0108 (2001) 014 (hep- 
ph/0106258), and hep-ph/011115(\ 

[7] A. Bandyopadhyay et al., Phys. Lett. B519 (2001) 83. ( hep-ph/0106264 ). 

[8] P. Creminelli, G. Signorelli and A. Strumia, hep-ph/01022 3^ version 2 from 10 July, 2001. 


[9] M.V. Garzelli and C. Giunti, ^ep-ph/0108191 


[10] P. 1. Krastev and A. Y. Smirnov, fiep-ph/0108177 

]11] M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B466 (1999) 415 ( hep-ex/990703'^ . 

]12] F. Boehm et al, Phys. Rev. Lett. 84 (2000) 3764, Phys. Rev. D62 (2000) 072002 and Phys. Rev. 
D64 (2001) 112001. 

]13] MINOS Collaboration, P. Adamson et al., NuMl-L-476 (March 1999); Web page: 
http://www. hep. anl.gov/ndk/hypertext/numi.htnil 

[14] A. Piepke [KamLAND Collaboration], Nucl. Phys. Proc. Suppl. 91 (2001) 99. 

[15] S.M. Bilenky et al, Phys. Lett. B 356 (1995) 273 and Phys. Rev. D 54 (1996) 1881. 

[16] S. M. Bilenky, C. Giunti and W. Grimus, Prog. Part. Nucl Phys. 43 (1999) 1 ( hep-ph/9812360^ . 

[17] B. Pontecorvo, Zh. Eksp. Teor. Fiz. 33 (1957) 549, and 34 (1958) 247. 

[18] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[19] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59 (1987) 671. 


[20] S. M. Bilenky, S. Pascoli and S. T. Petcov, Phys. Rev. D 64 (2001) 053010 (hep-ph/0102265). 

[21] P. Lipari, Phys. Rev. D61, 113004 (2000). 

[22] V. Barger et al, Phys. Rev. D62 (2000) 013004. 

[23] M. Freund et al, Nucl. Phys. B578 (2000) 27. 

[24] M.C. Banuls, G. Barenboim and J. Bernabeu, Phys. Lett. B531 (2001) 391; J. Bernabeu et al. 


hep-ph/0110071 


[25] R. Barbieri and A. Strumia, fiep-ph/001130 / . 

[26] S. Schoenert, Talk given at the TAUP2001 International Workshop, September 8 -12, 2001, Lab- 
oratori Nazionali del Gran Sasso, Italy. 


[27] A. Strumia and F. Vissani, hep-ph/0109172 


[28] S.M. Bilenky et al, Phys. Rev. D 54 (1996) 4432 (hep-ph/9604364). 

[29] S. M. Bilenky, D. Nicolo and S. T. Petcov, hep-ph/01 12216^ 

[30] G. L. Fogli. E. Lisi and A. Palazzo, hep-ph/01 0508(\ 


[31] M. C. Gonzalez-Garcia et al, Phys. Rev. D 63 (2001) 033005 ( jiep-ph/000935(/\) . 

[32] B. Pontecorvo, Chalk River Lab. Report PD-205 (1946). 


12 




































[33] For reviews see: G. Gratta, Nucl. Phys. Proc. Suppl. 85, 72 (2000) f \iep-ex/9905011 ), and 
C. Bemporad, G. Gratta and P. Vogel, fiep-ph/0107277 

[34[VI. Kopeikin, fiep-ph/01 1003(\ 

[35[ VI. Kopeikin, L. A. Mikaelyan and V V Sinev, Phys. Atom. Nucl. 64, 849 (2001) [Yad. Fiz. 64, 
914 (2001)[ ( fiep-ph/01102^ . 

[36[ See, e.g., C. Bemporad, G. Gratta and P. Vogel in 

[37[ M. C. Gonzalez-Garcia and C. Pena-Garay, hep-ph/0111432\ and references therein. 

[38[ V. Barger et ah, hep-ph/0110393. 

[39[ W. Pauli, Letter to Participants of a Physics Meeting in Tubingen (H. Geiger, L. Meitner, etc.), 
Germany, December 4, 1930. 


13 

















